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FOREWORD 


In  evaluating  any  marine  system  (ASW,  tracking,  .etc.), 
it  is  desirable  to  correlate  system  performance  with  en- 
vironiiiental  conditions.  Sound  velocity  is  often  the  only 
oceanographic  parameter  monitored  during  system  tests . 

This  report  describes  a  technique  for  determining  gross 
distributions  of  teii5)erature  and  salinity  from  sound  ve¬ 
locity  profiles .  This  work  was  acco^lished  under  Bureau 
of  Naval  Weapons  WEPTASK  Assignment  No.  RUTO-3E-OOO/219 
l/SP-099-03-02  and  RU22-2E-000/219  1/R004-03-01. 
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SUMMARY 


This  report  presents  ^technlqiie^for  the  interpretation 
of  sound  velocity  profiles^  in  terms  of  “basic  oceanographic 
parameters  of  tenperature ,  salinity  and  depth.  An  enpirical 
relationship  between  these  parameters  and  the  speed  of 
sound  in  sea  water  has  been  developed  by  Wayne  D.  .Wilson  of 
the  U.  S..  . Naval  Ordnance  Laboratory.  Using. Wilson's  equation 
as  a  tool,  graphs  have -been  prepared  which  show  the  rate  of 
change  of  sound  velocity  with  respect  to  temperature ,  .salin¬ 
ity,  pressure,  .or  any  combination  of  these  parameters.  Typ¬ 
ical  sound  velocity  profiles  are  also  presented  and  inter¬ 
preted  in  terms  of  the  environmental  conditions  which  pro¬ 
duce  them.  These  graphs  and  analyses  should  help  engineers 
to  infer  general  oceanographic  conditions  from  given,  sound 
velocity  profiles,  and  to  determine  the  strength  and  extent 
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nraRODUCTION 


Sound  velocity  Ig  not  the  only  oceanographic  parameter  which  affects 
the  performance  of  marine  systems.  Temperature  and  salinity  distributions 
are  often  as  useful  and  sometimes  more  important.  But,  after  conpleting 
a  series  of  acoustic  tests  in  which  only  sound  velocities  were  taken, 
test  engineers  -often,  discover  that  it  would  be  useful  to  have  some  idea 
of  the  gross  distributions  of  temperature  and  salinity  as  well.  This  re¬ 
port  presents  a  technigiie  for  inferring  the  gross  distributions  of  tenper- 
ature  and  salinity  from  sound  velocity  data  and  from  a  general  knowledge 
of  oceanographic  distributions. 

Wayne  D.  Wilson  of  the  U.  S.  Naval  Ordnance  Laboratory  has  developed 
an  equation  which  relates  tenperature ,  salinity  and  pressure  of  a  sea 
water  sample  with  the  speed  of  sound  propagation  through  the  sample.  ,A 
full  explanation  of  Wilson's  equation  in  terms  of  basic  physical  relation¬ 
ships  and  experimental  results  is  contained  in  NAVORD  Report  No.  6986 
(see  bibliography).  .For  purposes  of  this  report  Wilson's  equation  serves 
6is  a  tool  in  constructing  graphs  ■vdiich  can  be  used  by  engineers  to  calcu¬ 
late  variations  in  soimd  velocity  with  tenperature,  salinity  and  pressure, 
and  to  explain  typicsil  vertical  sound  velocity  profiles  in  terms  of  tem¬ 
perature  :6Uid  salinity  depth  distribution,. 


SOUND  VELOCITY  IN  SEA  WATER 


Wilson's  equation  for  the  speed  of  soiand  propagation  In  sea  water  is: 
V  =  li4l|-9.1^  +  A  Vrp  +  A  Vp  +  A  Vg  +  A  Vgijip. 

A  Vt  »  CiT  +  GgT^  +  G3T3  +  ChT^. 

AVp  -C5P  +  CgP^  +  C^P3  +  CqP^ 

A  Vg  -  09(3-35)  -+  Cio(S-35)^. 

A  Vgip  =  (S-35)[CiiT  +  +  C.13P  +  .Cx4P^  +  C15PT  +  CigPT^] 

+  PtCijT  +  +  Ci^T^  l  +  P^  [e2oT  +  C2iT^  ]  +  022^?^. 

In  this  equation:  :V  =  velocity  of  sound  in  meters  per  second 

T  m  tenperature  in  degrees  centigrade 
.,.;S  ■  salihity  in  grams  per  kilogram  of  sea  water 
/?  ■  pressure  lii  kilograms  pxer  square  centineter 
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•Cl  -  4.5721  p 
C2  -  -4.4531  X  10' J 

Co  -  -2. 6045  X  10;^ 

€4  -  T. 98.51  X  10"^^ 
C5  -  1.60272  X  10"-^ 
Cg  ■  1.0268  X  10*5 

Cj  -  3.5216  X  10-9  „ 

00  -  -3.3606  X  10"-^^ 
Cq  -  1., 39799 

CiQ  m  1.6920  X  10-^ 

■  -1.1244  X  10"2 


7.7711  X  lO-'i' 

7.7016  X  10-5 
1.2943  X  10-7 
3.1580  X  10-^ 
1.5790  X  10-9 

-1.8607  X  10"^ 
7.4812  X  10-9 
4.5283  X  10"° 

-2 . 5294  X  10-7 
1.8563  X  10-9 
-1.9646  X  10"^° 


The  speed  of  somd  propagation  in  sea  water  is  shown  to  vary  with  tem¬ 
perature,  .salinity,  pressure  and  a  combination  of  all  three.  To  determine 
these  variations,  the  first  step  is  to  obtain  the  partial  derivatives  of 
the  velocity  equation  with  respect  to  temperature,  .salinity  and  pressure. 

-  Cl  +  2C2T  +  3C3T2  +  4C4t3  +  Cii(S-35)  ■+  2Ci2(S-35.)T 

+  Ci5(S-35)P  +  2Ci6(S-35)PT  Ci^^P  +  2Ci8K?  +  3C19PT2 

+  C2/  +  2C2iP^T  +  C22^^ . 

11  -  C9  +  2Cio(S-35)  +  CiiT  +  CigT^  +  C13P  +  Ci4p2 

+  C15PT  +  CigPT^. 

-|1  -  C5  +  2C6P  +  30^?^  +  4C8P3  +  013(8-35)  +  2Ci4(S-35)P 

+  015(8-35  )T  +  Ci5(8-35)T^  +  Ci.j.T  +  CisT^  +  Ci^T^  +  2C2(fT 
+  2C2iPT^  +  3022?^^'. 

Prom  these  partial  derivatives  it  is  evident  that  the  variations  of 
velocity  wlth-tem^^rat'Uire^  salinity  and  pre^s'tod  are'-theinseiVeS  func-  " 
tionslb^  the...  teinperature ^  salinity  and  pressure.  As  there  is  no  convenient 
way  in  which  to  present  the  variation  of  a  quantity  with  respect  to  three 
variables;  the  variation  of  each  partial  de.rivatlve  with  respect  to  salinity 
was  investigated-.  (This  variation  is  represented  as  the  partial  derivative 
with  respect  to  seiinlty  of  the  temperature;  aalinity  and  presstire  partial 
derivatives  already  obtained. )  The  resulting  equations  are: 


(M) 


Cii  +  ^^12^  ^  ^15^  ^ 


2 
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1 . 

^Cbs/^ 

L 

as 

1  4 

i-i 

1  : 

-m. 

as 

C13  +  2C14P  +  C3_5T  -f-  C-^^1^. 


These  equations  indicate  that  the  rate  of  change  with  respect  to  salinity 
of  the  rates  of  change  of  sound  velocity  (with  respect  to  temperature,  salinity 
and  pressure)  is  a  function  of  temperature  and  pressure  only.  Using  this  re¬ 
lationship,  the  rate  of  change  of  sound  velocity  with  respect  to  tenperatiure 
(which,  as  already  noted,  is  a  function  of  ten5>erature ,  salinity  and  pressure) 
may  he  expressed  as  a  function  of  tenperat\jre ,  pressure  and  a  constant  refer¬ 
ence  salinity  (Sq)  plus  the  difference  between  the  actml  salinity  and  the 
reference  salinity  multiplied  by  the  rate  of  change  with  respect  to  salinity 
of  the  rate  of  change  of  sotmd  velocity  with  respect  to  temperatijre;.  The 
same  procedijre  may  be  used  to  express  the  rate  of  change  of  sound  velocity 
with  respect  to  salinity  or  pressure.  The  three  equations  thus  obtained  are: 


11  =Jl(T,S,P)  i(T,So,P)  +  (S-So)  — • 


-||  -52(T,S,P)  -f2(T,So,P)  +  (S-Sq) 


^  -  53(T,S,P)  -  f3(T,So,P)  +  (S-Sq) 


M 

as 

■ 


These  substitutions  are  j\astif ied  over  the  entire  range  of  salinity 
because  the  rates  of  change  of  sound  velocity  with  respect  to  temperature, 
salinity  and  pressure  are  all  linear  in  salinity.  The  reference  salinity 
selected  was  35  grams  per  kilogram  of  sea  water  -  a  common  salinity  in  the 
world's  oceans. 


VARIATION  WITH  TEMPHIATURE  ,  SALINITY  AND  PRESSURE* 


Figures  1  through  6  illustrate  graphically  the  general  characteristics 
of  variations  in  sound  velocity  with  respect  to  temperature ,  salinity  and 
pressure,  and  the  interaction  of  thfese  parameters.  Figvire  1  shows  the  rate 
of  chan^  of  sound  velocity  with  respect  to  temperatinre  (  3  V/dT)  as  a 
function  of  teii:q)erature  and  pressxire  (salinity  constant  at  35  o/pb).  It  is 
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evident  that  teii5>erature  has  a  far  greater  infl-uence  than  pressure  on. 
aV/8T.  .Moreover,  the  effect  of .  temperature  is  inverse;  S.V/aT  de¬ 
creases  with  increasing  temperatiire.  The  effect  of  pressure  on  S.V/^T 
appears  to  depend  on  the  temperature  .  At  0“C  3  V/8  T  decreases  with  in¬ 
creasing  pressure,  whereas  at  30°C  av/^T  increases  with  increasing  pres¬ 
sure.  The  effect  of  salinity  on  aV/aT  is  shown  in  f:lgure  2  :«is  a  function 
of  temperature  and  pressure.  8  (d  V/8  T)/8  S  varies  inversely  with  both 
te^5)erat^ire  and  pressure,  the  effect  of  pressxa^  being  greater  at  higher 
temperatures. 

Figure  3  shows  the  rate  of  change  of  sound  velocity  with  respect  to 
salinity  (SV/aS)  as  a  function  of  tenperature  and  pressure  (salinity 
again  constant  at  35  b/od).  3  V/  8  S  varies  inversely  with  temperature:. 

With  respect  to  pressure  a  V/ 3  S  is  almost  constant ,  ,althou^  there  is 
an  indication  that  it  builds  to  a  maximum  with  increasing  pressure  :and 
then  decreeises  as  pressure  continues  to  increase.  Ihe  pressure  at  which 
this  maximum  occurs  appears  to  vary  slightly, with  temperature.  .With  respect 
to  salinity  3V/3S  increases  linearly  at  the  rate  of  3.38^  x  10“'8  meters/ 
second  per  gram/kilogram  per  gram/kilogram  of  salinity. 

Figure  k  shows  the  rate  of  change  of  sound  velocity  with  respect  to 
pressure  (8  V/8P)  as  a  function  of  temperature  and  pressiu’e  (salinity 
again  constant  at  35  d/dd).  It  is  very  difficult  to  describe  the  vari¬ 
ations  shown  in  this  figure;  therefore,  the  same  material  is  presented 
in  another  form  in  figure  5.  In  this  form  3  V/3P  increases  with  Increasing 
pressure  for  temperat\ires  less  than  20“G.  The  rate  of  increase,  however, 
decreases  with  Increasing  temperatitre .  For  temperatures  above  2,0“e  and 
pressures  greater  than  200  kilograms  per  sq,uare  centimeter,  the  lines  of 
constant  3  V’/3  P  reverse  their  curvature  and  become  concave  upward.  No 
explanation  of  this  will  be  offered  because  the  purpose  of  this  report  is 
to  explain  the  environment,  not  the  physics  behind  Wilson's  eq.xiatlon. 

.Also,  this  area  (>20®C  and  100  k^i^cmr)  is  outside  the  limits  of  natural 
variability  in  the  world's  oceans.  The  effect  of  salinity, on  3,V/3P, 
shown  in  figure  6  as  a  function  of  temperatxu'e  and  pressure,  is  very  small. 

.  .From  these  graphs  (figures  1  throu^^  6,)  the  rate  of  change  of  /sound 
velocity  with  respect  to  temperature,  salinity,  .or  pressure  may  be  obtained 
(as  follows )  for  any  comiblnation  of  temperature ,  .salinity  .and  pressiire  ob¬ 
served  in  the  ocean  to  a  depth  of  5000  meters .  ' 

,1..  Locate  the  graph  showing  the  rate  of  change  of  sound  velocity 
•with  respect  to  the  desired  parameter  .(figure  1,  >3,  h-  or  5)* 

2.  Locate  the  corresponding  graph  for  salinity  corrections  (figure 
2  or  6).  The  salinity  correction  to  the  rate  of  change  of  somd  velocity 
with  respect  to  salinity  (3.V/33)  is  etpal  to  3*384  x  10"3  meters/secdnd 
per  gram/kilogram  of  selinity  per  grami/kilogram  of  salinity  at  all  temper¬ 
atures  and  salinity. 
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3:.  Read  the  value’  for  rate  of  chEuige  of  sound  velocity  at  the 
desired  tei!$erature  and  pressure  from  the  graph  selected  in  step  1. 

4.  Determine  the  salinity  correction  factor  at  the  desired  tem¬ 
perature  ;ahd  salinity  from  the  graph  selected  in  step  2. 

5.  Subtract  the  reference  salinity  (35  o/oo)  from  the  desired 
salinity  and  multiply  this  difference  by  the  salinity  correction  factor 
determined  in  step  4. 

6.  Add  the  salinity  correction  determined  in  step  5  to  the  value 
for  rate  of  change  of  sound  velocity  obtained  in  step  3. 


INTERPRETING  SOUND  VELOCITY  PROFILES 


Standard  vertical  sound  velocity  profiles  express  soxuid  velocity  as  a 
function  of  depth.  The  rate  of  cheuage  of  sound  velocity  with  respect  to 
depth  ( a  v/3^  may  be  considered  as  the  sum  of  the  rates  of  change  of 
sound  velocity  with  respect  to  temperature,  salinity  and  pressure,  each 
multiplied  by  the  rate  of  change  of  the  same  parameter  with  respect  to  depth. 


The  partial  derivative  notation  has  been  \ised  to  indicate  that  temperatiare , 
salinity,  pressure  and  sound  velocity  may  vary  in  time  or  spewe  at  constant 
depth. 

In  considering  sound  velocity  profiles  it  should  be  noted  that  dT/jz; 
and  3S/dZ  vary  over  a  wide  rar.ge,  whereas  the  variations  in  8  P/a  2 
are  quite  small  so  that  in  nany  cases  a  constant  value  may  be  used.  Appendix 
A  contains  a  coiqilete  disc\ission  of  the  factors  that  contribute  to  pressure 
at  any  depth.  For  general  interpretation  of  sound  velocity  profiles,  how¬ 
ever,  it  is  not  necessary  to  consider  variations  of  sea  water  composition, 
con5)ressibility  of  sea  water,  variation  of  gravitational  forces  with  depth, 
or  other  such  factors.  ap/ai  may  therefore  be  approximated  as  0.10 
Kgm/cm^  per  meter  of  depth. 

The  sound  velocity  profiles  in  figures  f  and  8  are  indicative  of  the 
effect  of  pressure  variation  with  depth  on  sound  velocity  for  various  tem¬ 
peratures  and  salinities.  The  displacement  of  the  profiles  in  figure  7  is 
due  to  differences  in  temperature  (salinity  constant  at  35  b/oo).  Although 
the  four  profiles  appear  similar,  they  vary  sli^tly  due  to  differences  in 
8V/8P  at  different  temperatures.  The  profiles  are  really  curves  and  not 
straight  lines  because  at  constant  temperature  av/a  P  varies  with  pres sixre. 


I 
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The  variations  of  av/a|P  with  teiiQ)erature  and  pressure  are,  however,  quite 
small;  hence  the  curves  appear  to  have  the  same  constant  slope.  Likewise,  w 

in  figure  8  the  profiles  for  all  three  values  of  salinity  (temperature  con¬ 
stant  at  20“C)  appear  to  ha^'/e  similar  and  constant  slopes.  There  are 
actually  small  differences  due  to  variations  in  9  V/  ^  P  with  salinity 
( teniperature  and  pressure  constant)  and  with  pressure  (temperature  and 
salinity  constant).  The  profiles  in  figures  T.and  8  show  that  even  if 
tenqjerature  and  salinity  are  consteuit,  sound  velocity  increases  by  .approx¬ 
imately  17  neters  per  second 'for  every  1000  meters  of  depth. 

The  slope  of  the  sound  velocity  profile  will  vary  with  salinity  .and 
tenperature  gradients.  Salinity  gradients  are  considered  separately .from 
temperature  gradients  to  illustrate  the  effect  of  each  independently  froia 
the  other . 

Figure  9  shows  the  sound  velocity  profiles  for  three  constant  salinity 
gradients,  starting  with  a  salinity  of  35  o/oo  on  the  surface  (temperature 
constant  at  20‘’C).  The  profiles  show  that  a  salinity  gradient  stronger  than 
-.01  o/oo  per  meter  is  required  to  offset  the  increase  in  sound  velocity  d\ie 
to  pressure  (depth).  The  strength  of  the  salinity  gradient  required  for 
isovelocity  (constant  velocity)  may  be  calculated  from  the  g  v/ Q  li  equation: 


For  isovelocity  conditions :  9V/9Z  *»  0. 

For  constant  temperature  :  3T/9Z  =  0. 

Therefore.:  (3V/3T)  (9T/9^  )  «  0. 

Making  these  substitutions : 


Referring  to  figure  3:  for  a  salinity  of  35  o/oo  and  a  temperature  of  20“C', 
a  V/9  S  c;  1.17  meters/second  per  gram/kilogram  of  sedinity  (for  pressures 
less  than  100  kgm/cm?).  Referring  to  figure  5:  for  the  same  conditions, 
ay/3P  —  0.l60  neters/second  per  gram/kilogram  (for  pressures  less  than 
100  kgm/cn^).  .Substituting  these  values  and  the  approximation  3  P/9Z  ■ 
0.. 10  into  the  equation: 


0 


1.17  +  (•16q)  (.10). 
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This  yields;  ^:S/5E=  0.0137. 

'Eierefore^  tinder  the  conditions  given>  a  salinity  gradient  stronger  than 
-0.013T  o/oo  per  meter  is  required  to  cause  sound  velocity, to  decrease  with 
depth. 


The  effect  of  constant  temperature  gradients  on  the  sound  velocity  pro¬ 
file  is  shown  in  figure  10  (salinity  constant  at  35  o/oo).  .For  a  temperature 
of  20®C  at  the  surface ^  .a  ten^rature  gradient  of  -0.00-5:°C  per  meter  will 
cause  approximate  isovelocity.  The  fact  that  sound  velocity  decreases  with 
increasing  depth  even  though  the  teii5>erature  gradient  remains  constant  is 
due  to  an  increase  in  3  V/ST  with  decreasing  temperature,  causing  the 
(^V/3T)  (  a  t/ 5  Z  )  term  to  outwei^  the  (5V/5P)  (aP/aZ  )  term, 
which  is  not  growing  so  rapidly. 

Thus  far,  variations  in  sound  velocity  profiles  have  been  .considered 
in  terms  of  constant  teiip)erature  or  salinity  gradients .  By  considering 
the  sum  of  one  or  more  temperatinre  and  salinity  gradients,  it  is  possible 
to  form  every, type  of  sound  velocity  profile  found  in  the-  ocean.  There 
are  several  general  types  of  profiles  which,  with  small  variations,  .account 
for  most  conditions  found.  Each  type  of  profile  must  be  considered,  in  de¬ 
tail  if  general  conclusions  eire  to  be  drawn  concerning  the  relative  effects 
of  temperature,  salinity, and  pressure  on  sound  velocity. 

Single  Minimum 

The  most  frequently  observed  type  of  sound  velocity  profile  is  the  single 
minimum,  in  which  a  single  minimum  sound  velocity  occurs  at  some  depth  (figure 
ll);.  This  type  of  profile  may  be  explained  in  terms  of  a  very  simple  set  of 
oceanographic  conditions.  If  salinity  is  assumed  to  be  constant,  a  single 
minimum  profile  can  arise  from  a  ixniform  temperature  gradient,  extending 
from,  the  surface  to  the  depth  at  which  minimum  sdxnd  velocity  occurs  and 
terminating  at  that  depth  in  a  consttot  temperatvife  layer.  .A  check  of  the 
curves  and  gradients  already  presented  will  verify  . this  intei^retation  and 
allow  calculation  of  the  gradients  involved.  This  analysis  presents  the 
overall  picture  without  taking  into  account  very  slight  salinity  gradients 
or  slight  changes  in  the  temperature  gradient. 

The  single  minimum  profile  is  found  in  most  ocean  areas  because  a  tem¬ 
perature  gradient  from  the  surface  to  some  depth  does  exist.  This  gradient 
is  d\ae  to  deep  water,  formed  in  the  arctic  or  antartic,  which  Is  colder  than 
the  water  found  in  miost  areas  of  the  ocean.  Although  salinity  is  not  con- 
-stant  in  open  ocean  areas ,  the  range  of  variation  is  small  and  sallhlty 
gradients  vary  but  slightly.  ^ 
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Surface  Channel 


III  all  sound  velocity  profiles  (noticeable  or  not)  there  is  a  surface 
sound,  channel  -  -aii  area  of  increasing  sound  velocity  (figure  12)  --which 
•results  fromi  isothermal,  and  isohELllne  (constant  ten^^rature  and  salinity) 
conditions  at  the  surface:.  .This  surface  layer  of  isothermal  euid  isohaline 
water  is  present  at  all  times  in  all  ocean  areas.  Thickness  of  the  layer 
varies  from  a  few  millimeters  to  several  hundred  meters .  These  conditions 
result  from  the:  dynamic  mixing  of  water  and  the  layer  is  often  referred  to 
as  the  mixed  layer.  , Tenqperature  and  salinity  gradients  for  which  (3V/9  T) 
(^T/3^  )  +  (av/3S)  Os/32  )  <  (av/ap)  (a;p/az  )  may  .also  con- 
-tribute  to  the  thickness  of  this, •s\arf ace  sound  channel. 

Surf ace , Channel  and  Single  Minimum 

The  surface  sound  channel  exists  on  all  sound  velocity  profiles.  .Figure 
13  shows  the  surface  soimd  charmel  on  a  single  minimum  profile.  .The  temper - 
-atiire  gradient  is  in  three  simple  sections:  a  constant  temperat-uTe  layer 
-at  the  siarf.ace,  ,a  thermocline  or  layer  in  which  temperature  decreases  with 
depth,  .and  a  constsuit  -tenperature  layer  at  the  bottom.  Salinity  is  consid¬ 
ered  as  constant  and  has  no  effect  on  the  profile.  .This  is  a  very  common 
■type  of  profile  becaiase  the  conditions  described  are  common  to  all  water 
masses. 


Double  .Minimum 

Salinity  is  not  always  constant.  In  aresis  near  -shore,  in  bays,  .and 
near  large  ice  masses,  .salinity  may  vary  considerably  with  depth.  .When  -a 
large  positi-ye  change  in  salinity  occurs  in  the  thermocline,  the  sound  veloc¬ 
ity  profile  may  have,  two  minimums  (figures  l4  and  15 ).  Sound  -velocity,  -which 
■is  decreasing  with  depth  in  the  thermocline  ,  .increases  with  depth  in  the 
halocline  (salinity  gradient).  .Af-ter  passing  throu^  the  halocline,  the 
sound- velocity  again  decreases  with  depth  lantil  the  thermocline  -teiminates , 
,at  which  point  it  increases  dvie  to  the  pressure  effect-. 


There  are  two  types  of  profiles  with  a  double  minimum:  the  shallow 
minimum  and  thb  deep  minimum.  .The.  conditions  for  both  are  the  same  -  -a 
layer  -of  relati-vely,  fresh  wa-ber  at  the  s\jrface:.  The  difference  is  due  to 
the  volume  of  fresh  -water  or  thickness  of  the  fresh  -water  layer A  large 
vol-ume  of  fresh  water  -  e-g..  ,  near  a  melting  ice  pack  -  produces  a  thick 
stirface  layer,  .and  the  halocline  occurs  deep  in  the  thermocline  (figure  l4). 
Under  these  conditions  the  shallow  minimum  has  a  lesser  -value  (in  terms  of 
:sound  -velocity)  than  the  deep  minimum,,  .because  there  is  insufficient  thermo- 
cllne  left  below  the  halocline.  to  produce  more  of  a  decrease  in  sound  velocity 
than  has  already  ocC-urred  at  the  sheLllow  minimum-.  On  the  other  haai^  a  small 
volume  of  fresh  water  -  -e.g. ,  near  shore  areas  with  fresh  -wa-ter  run-off  - 
produces  a  thin  svcrface  layer  so  that  the  halocline  occ\irs  :Shallow  in  the 
thermocline.  Under  these  conditions  the  deep  mihimum  has  -a  lesser  -value 
(figure  15)i 
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Kink 

If  a  hi^.  saline  layer  overlays  a  lass  saline  layer,,  the  result  is  a 
kink  type  of  sound  velocity  profile  (figure  16).  Such  .conditions  may  be 
fo\md  in  tropical  areas  where  evaporation  exceeds  precipitation,  or  near 
a  river  rvin-off  during,  winter  where  there  is  an  intrusion  of  cold  fresh 
water  which  is  denser  than  the  surface  waters.  This  decrease  in  salinity 
results  in  an  additional  decrease  in  sound  velocity  in  the  halocline, 
giving  the  profile  a  kinked  appearance-. 


GONCLUSION 


This  report  has  analjraed  the  effect  on  sound  velocity  of  changes  in 
temperature,  salinity  and  pressure,  based  on  Wilson's  equation  for  the 
variation  of  sound  velocity  in  sea  water.  To  illustrate  this  variation  a 
series  of  graphs  have  been  prepared  which  show  the  rate  of  change  of  soimd 
velocity  with  respect  to  ten^ierature ,  salinity,  pressure,  or  any  combination 
of  these  parameters ,  Typical  sound  velocity  profiles  have  been  presented 
and  interpreted  in  terms  of  the  environmental  conditions  that  produce  them. 

By  using  these  graphs  and  typical  profiles,  -engineers  can  infer  the 
general  oceanographic  conditions  from  sound  velocity  profiles  and  deter¬ 
mine  the  strength  and  extent  of  teiip)erat\u'e  and  salinity  gradients. 
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Figure  1  •  Rate  of  Change  of  Sound  Velocity  with  Respect  to  temperature. 


01125 


Figure  2.  Rate  of  Change  of  Sound  Velocity  with  Respect  to  Temperature:  Salinity  Correction. 
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FIGURE  3 


Figure  3*  Rate  of  Change  of  Sound  Velocity  with  Respect  to  Salinity. 


Figure  4.  Rate  of  Change  of  Sound  Velocity  wjth  Respect  to  Pressure. 


RATE  OF  CHANGE  OF  VELOCITY  WITH  PRESSURE  FOR  SALINITY 


Rate  of  Chanae  of  Sound  Velocity  with  Respect  to  Pressure. 
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Figure  6.  Rate  of  Change  of  Sound  Velocity  with  Respect  to  Pressure:  Salinity  Correction* 


SOUND  VELOCITY  PROFILES  AT  SALINITY  OF  35*/. 


Fi^re  7.  Sound  Velocity  Profiles  for  Various  Constant  Temperatures 
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Figure  8.  Sound  Velocity  Profiles  for  Various  Constant  Salinities 


Sound  Velocity  Profiles  for  Various  Salinity  Gradients. 
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Ftgur*  11.  Sound  Volocity  Profile: 
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Rguro  12.  Sound  Velocity  Profile:  Surface  Channel. 
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Figure  13.  Sound  Velocity  Profile:  Surface  Channel  and  Single  Mirtimum. 
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Fijpire  14.  Sound  Velocity  Profile:  Double  Minimum  Shallow. 
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.APPENDIX  A 


(Material  in  this  appendix  is  based  on  NUGS  Technical  Memorandian 
No.  2&h,  ."Processing  Gceanogr'apbic  Station  Data  Using, Rapid  Data-Handling 
Eq.uipiiient " ^  C.  .F.  Morey,,  U.  S.  Naval  Underwater  Ordnance  Station,  .August  I962. 
For  references  to  formulas  and  publications  see  this  earlier  report-. ) 

To  coH5)'ate.  the  velocity  of  sound  in  sea  vater  according  to  the-  formulas 
developed  by  Wayne  D.  .Wilson  (which  give  velocity  as  a  function  of  ten^ier- 
ature,  .salinity .and  pressure),  the  pressure  must  be  calculated  from  the 
observed  oceanographic  station  data.  The  pressure  at  any  point  in  a  fluid 
is  the  sum  of  the  atmospheric  arid  hydrostatic- pressures-.  Hydrostatic  pres- 
■sijre  is  given  by: 

A? 

Pn  ^g  d'S 

where :  •Pn  o  hydrostaMc  pressure  in  .  ,newtona/ciri^, 

»  density  in  kgm/cm’^. 

g  a  gravltatioheU.  acceleration  in  cm/sec^. 

2  depth  in  nKsters . 

A  method  for  evaluating  this  equation  for  the  pressure  in  a  non-homogeneous 
compressible  fluid  in  a  non-uniform  gravitational,  field  haa  been  e:q)ressed  by 
N.  -P.  .Fofonoff  and  N.  A*  Pedersen.  This  method  has  been  adopted  vith  modi¬ 
fications  by, the  Naval  Underwater  Ordnance  Station  for  calculating  press\ire 
in  the  ocean  from  oceanographic  station  data-. 

The  geographic  variation  of  gravitational  acceleration  is  best  expressed 
by  the  formula  accepted,  by  .  the  General  Assembly  of  the  International,  Union 
of  Geodesy. and  Geophysics  at  Stockholm^  Sweden  in  1930«  This  formula  is: 

go  =  978.049  (1  +  0-. 005288.4  sin^j)  -  -0.0000059  s.in^(|)) 
where:  go  «  gravitational  acceleration  at  sea  level  in  cm/sec^.. 

|)  •«  latitude  in,  degrees. 

The  variation  of  gravitational  ;accel6ritionitfithodepth.lias^been  Igiven^by 
Sverdeep  as :  * 

g  ,=-  .gQ..+- 2':.,202 '  x  io'-^2 

This  expression  has  been  verified  by  me.aSureinenbs  made  in  the  hathyscaph 
TRIESTE  by  K.  Mackenzie. 
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V.  .W.  Ektnan  has  expressed  the  variation  in  density  as  a  function  of 
temperature,  .salinity  and  pressure.  Hie  complete  functional  relationship 
hetwen  these  quantities,  .which  is  very  con^lex,  is  eiqpressed  here  in 
terms  of  functions  which  are  frequently  computed  in  oceeuiographic  station 
work. 

<^S,T,P 

vhere:  /^s,T,P  ®  density  as  a  function  of  salinity,  ten?)erature  and  pressure. 

T  P  “  specific  volume  as  a  function  of  salinity,  teii5)erature 
and  pressure . 

The- specif  ic  volume  is  expressed  eis  : 


- 


CCq  m  -p.m  ttg  „  Q  _  P  QU  „  Q  10~9r^.‘P§.^.  ■?,  ^  +  (A  +  B) 

S,T,P  S,,T,0  -^,T,0  ^  .3_.g2  - 

where:  A  =  P(l0"^)[l05.5  +  9.5T'  -  0.158t2]  -  1.5p2t(10"®) 

-  ,[227  +  28.33T  -  0.551T^  +  O.OO^T^]. 

B  -  -  ^  [147.3  -  -2.72T  +  0.04T^  ^  P(10"^)(32:.4  -  O.87T  +  0.02T^)] 

+  (££^iJS)  [4.5  .  O.IT  -  P(l0-^)(l.8  -  0..,o6t)] 

K  =  -28.134 


a. 


1  +  o^(io-3) 

0^  »  Zip  +  (oq  +  0.1324) [1  -  Mp  +  NjCOq  -  0.1324)] 


T  -h  28  3 
T  +  67.26 J 

=  T( 4.7867  -  O.O98185T  +  0.0010843T^)10"^ 

%  T(18.03  -  o.8i64t  +  0. . 016671^ )io‘^ 

Pq  b  -0.06,9  +  i.47o8ci  -  q.ooi^tg-*^  +  o.oooo3§>8c-e^ 


Erp 


(T  -  3. 
503.570 


4, 
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The  expression  for  calculating  the  density,  which  ,is  req.uired  for  the 
pressure  caiculation,  .itself  req;Uires  that  pressure  he  known.  This  dif¬ 
ficulty  may  he  overcome  hy  an  iterative  process  commonly  used  in  ma,the=- 
matics.  An  approximation  ‘ of  the  pressure,  is  made  in  order  to  determine 
density  and  then  pressure.  The  computed  pressure  is  compared  with  the  ap¬ 
proximation  used.  If  they  do  not  agree  within  certain  limits,  the  process 
is  repeated  using  the  confuted  pressure  as  the  new  approximation  to  deter¬ 
mine  again  the  density  and  pressure.  .This  proce.ss  is  repeated  until  agree¬ 
ment  is  reached  between  the  approximated  and  computed  pressures. 

The  first  approximation  for  pressure  is  computed, -.as: 

P  =  r  g,(l  +‘^11,  X  10"'3  )d2 
J  o 

vhere?.:  »  average  aj. 

Conrputatlon  of  the  pressure  is  further  coir5)licated  hy.  the  fact  that 
oceanographic  station  data  consist  of  a  number  of  unique  saoiples  taken  in 
the  vater  column,  .and  not  a  continuous  record  of  temperature,  .salinity  and 
depth.  With  this  data  the  trapezoidal  rule  is  used  successively  for  evalu- 
■ating  the  pressure  integral,; 


Using,  this  process  the  pressure  !  approximation  becomes  : 


Pi 


»  Pi-1  + 


i/kjr2  +  (< 


)10.^ 

■  ^i^l 


-3 


(gi 


.+  g 


i-1 


)( Z,  -  z 


Once  the  hydrostatic  pressure  has  been  determined  in  newtons/cm?,  the 
pressure  for  coii5)utation  of  sound  velocity  may  he  obtained  simply  hy  adding 
the  atmospheric  pressure  (l  millibar  *»  10"^  newtons/cm^ ) .  This  pressure 
must  he  converted  to  kilo^ams/cm^  hy  multiplying  by  1/9.80665  kilogrEp..  per 
newton.  Wilson's  formula  for  sound  velocity  can  then  he  -used: 


V  -  lW9a4  +  A  Vt  +  A  Vp  +  A  Vg  +  A  VsTP- 

A  Vqi  =  4.5T21T  -  4.4531  X  10"2t2  .  -2., 6045  X  10‘M  +  7*9851  X  10"^T^. 
A  Vp  »  1.60272  X  10“%  +  1.0268  X  10“%^  +  3.5216  x;  10“%3  .3.^3606 

X  10"^%^. 
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AVg  =.  1.39799(8-35)  +  1.6920  x  10-3(3-35)2. 

A  ^STP  -  (S"35)  (-  1.1241  X  10-2t  +  7.7711  x  +  7.70l6  x  10’^ 

■+  1.2943  X  10-'^p2  +  3.1580  X  10-®PT  +  1.5790  X  10-9pr2) 

+  P(-  1.8607  X  10"^:  +  7.4812  X  10"^t2  +  4.5283  X  10"8t3) 

+  p2(-  2.5294  X  lO-'^T  +  1.8563  X  10-9t2)  +  p3(.  1.9646  X  10"^%). 
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